Dynamical vanishing of the order parameter in a fermionic condensate 
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We analyze the dynamics of a condensate of ultra-cold atomic fermions following an abrupt change 
of the pairing strength. At long times, the system goes to a non- stationary steady state, which we 
determine exactly. The superfluid order parameter asymptotes to a constant value. We show that 
the order parameter vanishes when the pairing strength is decreased below a certain critical value. 
In this case, the steady state of the system combines properties of normal and superfluid states - 
the gap and the condensate fraction vanish, while the superfluid density is nonzero. 



Recently, several remarkable experiments have demon- 
strated Cooper pairing in cold atomic Fermi gasespJ-Q. 
Key signatures of a paired state - condensation of Cooper 
pairs0,0] and the pairing gap0 have been observed. In 
addition, trapped gases provide a unique tool to explore 
aspects of fermion pairing normally inaccessible in super- 
conductors. One of the most exciting prospects is a study 
of far from equilibrium coherent dynamics of fermionic 
condensates [5^3 ; made possible due to the precise exper- 
imental control over interactions between atoms [irlfiij|. 
The dynamics can be initiated by quickly changing the 
pairing strength with external magnetic field. 

In the present paper, we determine the time evolution 
of a fermionic condensate in response to a sudden change 
of interaction strength. Initially, the gas is in equilibrium 
at zero temperature on the BCS side of the Feshbach 
resonance with a coupling constant gi > 0. At t = the 
coupling is suddenly changed to a smaller value gf > 
on the same side of the resonance, gi — > gf, Fig. ^ (inset). 
Ground states of the system at the old, gi, and new, gf, 
values of the coupling are characterized by corresponding 
BCS gaps, A, and A / , respectively. We consider the case 
Ai > Af. It has been shown previously that following the 
change of coupling, the time-dependent order parameter 
A(t) asymptotes to a constant value |A(t)| — > Aoo on 
a timescale ta = l/A^. Here we evaluate in terms 
of Ai and Af. 

We show that when the coupling is decreased below 
a certain critical value, Aqq vanishes, Fig. ^ On a ta 
timescale the system goes to a steady non- stationary 
state that combines properties of normal and superfluid 
states in a peculiar way. For example, the gap vanishes, 
while the superfluid density remains finite. Provided the 
system is continuously cooled, the BCS ground state with 
a gap A^- is reached on the energy relaxation timescale 
t £ , which is typically much larger than ta- Experimen- 
tal signatures of the novel state include the absence of 
the gap in rf absorbtion spectrum and zero condensate 
fraction after a fast projection onto the Bose-Einstein 
Condensation (BEC) side (see below). 

At times t <C r £ , dynamics of the condensate in 
the weak coupling regime can be described by the 
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FIG. 1: Time evolution of the BCS order parameter A(i) 
following an abrupt change of the coupling constant. The 
coupling is changed by varying the magnetic field on the BCS 
side of the Feshbach resonance (inset). The order parameter 
exponentially decays to zero (solid line) at times t ~ 1/Ai 
and the system goes into a gapless steady state. If the sys- 
tem is continuously cooled, the order parameter recovers to 
its equilibrium value Af (dashed) at times longer than the 
quasiparticle relaxation time r e . 



BCS model. Here we are interested in the thermo- 
dynamic limit, in which case one can use the BCS 
mean- field approach 12J. Using Anderson's pseudospin 
representation j one can describe the mean-field evo- 
lution by a classical spin 
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where Cj are single-particle energies relative to the Fermi 
level and s^ 1 = ± is?. The summation in Eq. is 
over \ej\ < Ep, where Ep is the Fermi energy. Dynamical 
variables Sj are vectors of fixed length, |s 3 | = 1/2. The 
BCS order parameter is A(t) — A x — iA y = gJ2 j 
Equations of motion for classical spins s^- are 
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Components of spins are related to Bogoliubov ampli- 
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determine Aoo, we use the exact solution for the dy- 
namics of the BCS model[||. Consider the following 
vector function of an auxiliary parameter u: L(it) = 
(3) — z/g + J2j s j/( u ~ € j)> where z is a unit vector along 



At t = the system is in the ground state with gap 
Aj. The ground state is obtained by aligning each spin 
Sj antiparallel to its "magnetic" field bj jl2j in order to 
minimize the total energy Q for coupling constant g = gi 
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and s?(i = 0) = 0. 

At t = the coupling is changed, gi — > <?/, and the 
initial spin configuration @ is no longer an equilibrium 
for i > 0. To determine the time evolution of the system, 
one has to solve Eqs. @ with initial conditions (0J. 

We start with a linear analysis. Solving Eqs. J5J (with 
g = gj) linearized around the spin configuration J3J, we 
obtain up to terms of order SA/Af 
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where SA = A f - Aj, w(e) = + Aj, and h(e) = 

sinh _1 (e/A/). In Eq. JSJ), besides the continuum limit, 
we took the weak coupling limit Ep/Ai — > oo. 

The long time behavior of the order parameter in the 
linear approximation is obtained from Eq. © by station- 
ary phase method (see also Ref. 
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At times f > ta the gap approaches a constant value 
Aoo = A f to order SA/A f . 

Even though the gap is constant, the state of the sys- 
tem is non-stationary 8]. According to Eq. J5J), at large 
times each spin s^- = s(ej) precesses in its own con- 
stant field hj — (—2A oo ,0,2ej). For example, for the 
x-component of spins we derive from the linearized equa- 
tions of motion, 
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However, the gap A(t) = g^2j s j(t) contains oscillations 
with many different frequencies. At large times they go 
out of phase and cancel out in the continuum limit. 

In the non-linear case it can be shown that the gap de- 
cays to a constant A^ by a similar mechanism [ij. To 



the z-axis. Using Eq. (0, we obtain dL 2 (u)/dt = 0, i.e. 
L 2 (it) is conserved for any u. 

To determine the steady state gap A^, we evaluate 
L 2 (it) at t = for the initial spin configuration 10} and at 
t 3> Ta, when each spin s(e) precesses in a constant field. 
Matching the two expressions, we derive for Af < Aj 

Af = Aj exp(— retan/c/2), A^ = Aj cos k, (8) 

where < K < 7r/2. This equation determines Aoo in 
terms of Aj and Af, see Fig- El 

We make several observations. The steady state gap 
reaches its maximum Aoo = Af at Af/Ai = 1. Other- 
wise, Aoo < Af. Expanding around the maximum, we 
find Aoo = Af — (5A) 2 /6Af up to terms of higher order 
in SA/Af. The linear term vanishes in agreement with 
Eq. @. 




FIG. 2: After a sudden change of the BCS coupling constant, 
the order parameter A(t) (inset) saturates to a constant value 
Aoo. The plot shows the steady state gap in units of Aj 
as a function of the ratio A//Aj. The exact result given by 
Eq. (JSJ is compared to numerical solution of Eq. J5j. Note 
that Aoo = for A f < 0.2A,. 

Most interestingly, we see from Eq. (JSJ that the order 
parameter vanishes in the steady state, A m = 0, when 
K = tt/2, i.e. for (A//Aj) c = e^/ 2 w 0.21. If A//Aj 
is below this critical ratio, Eq. JSJ) has no solutions - the 
system goes into a gapless steady state, Aoo = 0, for 
Af < 0.21Aj. We also note from Eq. © and Fig.^that 
Aqo has a cusp at A//Aj = e _7r / 2 , i.e. the transition to 
the gapless state is second order. 

To gain further insight into properties of the gapless 
state, we need to know the spin configuration in this 
state. At large times spin s(e) rotates in a constant field 
b(e) = — 2AqoX + 2ez, where x and z are unit vectors 
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along the x and z axes, respectively. In the gapless case, 
b(e) = 2ez. The spin configuration can be characterized 
by the angle 9(e) between s(e) and — b(e). This angle 
can be determined by matching the conserved quantity 
L 2 (u) at t = and t > ta , 
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FIG. 3: Plot of the average occupation number (ra(e)) as a 
function of energy in the gapless state for various values of 
Af/Ai < e _7r,/2 . The exact result given by Eq. Q is com- 
pared to numerical solution of Eq. (|5J . Note that the distribu- 
tion function around the Fermi level is smeared over a width 
Se oc Aj. The inset shows a spin s(e+) at energy < e+ <C A; 
just above and below the critical point. 

Note from Eq. J2J) that in the steady state the compo- 
nent of spin s(e) along the field, sii(e) = — cos6*(e)/2, is 
constant, while the one perpendicular to the field, |sj_| = 
sin 9(e) /2 rotates with a frequency uo(e) = 2^/ e 2 + A 2 ^. 
The approach to the steady state can be studied by lin- 
earizing Eq. (J2J around this state. In particular, we ob- 
tain the asymptotic behavior of A(i) 
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Here a = — cos p and ir/2 < p < ir is the solution of 
p = ln(Ai/A/) cot(p/2). The parameter a has a prop- 
erty a — > when Af/Ai — * e~™/ 2 and a — > 1 when 



Af/Ai —> 0. The coefficient a in Eq. is time- 

independent, A(t) and B(t) in Eq. I|12|) are decaying 
power laws, A(t),B(t) oc l/t v with 1/2 < v < 2. 

We conclude that the decay law of A(t) changes from 
power law to exponential as we cross the critical point. 
Above the critical point spins rotate with frequencies 
uj(e) — 2yje 2 + A 2 ^. The inverse square root decay 
and oscillations with frequency 2Aoo in Eq. Ijll|l are 
due to the square root singularity in the spectral den- 
sity de/duj oc u) /^uj 2 - 4A 2 X) , cf. Eqs. ffi . Below the 
critical point Aqo vanishes, uj(e) — 2e, and the square 
root anomaly disappears. 

What happens to the spin configuration as Af/Ai is 
varied across the critical point? Consider a spin at energy 
e + <C Aj just above the Fermi energy. Using Eq. we 
obtain sin6»(e + ) w sin6>(0) = 21n(A l /A / )/7r for A f > 
e-^^Ai and sin0(O) = 1 otherwise. The field b(e+) = 
— 2A 00 x+2e_|_z is along the x axis above the critical point 
(Aqo ^> e + ) and along the z axis below (Aoo =0). At 
Ai — Af we have sin 0(0) = 0, i.e. the spin is parallel to 
the x axis. As Af/Ai decreases, the x component of the 
spin also decreases until it vanishes at the critical point 
Af = e _7r / 2 Aj. Just above it the spin lies in the yz plane 
and rotates around the x axis. Below the critical point 
it rotates in the xy plane around the z axis (Fig. |2J). 

In the gapless steady state each spin rotates around 
the z axis. Its z component is time-independent, s z (e) — 
— cos9(e)/2. It is related to the average occupation num- 
ber per fermion species (n(e)) at energy e as (n(e)) — 
s z (e) + l/2 [13, i.e. (n(e)) = (I- cos 9(e))/ 2. The distri- 
bution function (n(e)) is smeared near the Fermi energy 
over a width 5e oc Aj (Fig. [3J . Note that in this respect 
the gapless state is similar to the BCS ground state - 
the smearing over a width Se oc Ai due to interactions is 
also present in the ground state distribution (@J. In the 
normal state s z (e) — — sgn e/2 and (n(e)) = 0(— e). 

Given relation one can reconstruct the time- 

dependent condensate wave function, and evaluate nor- 
mal and anomalous correlation functions, e.g. 
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where c ea (c\ a ) annihilates (creates) a fermion of one of 
the two species a =|, { on energy level e. Note that even 
though Aqo vanishes, anomalous averages do not. 

Next, we determine the supcrfluid density n s in the 
gapless state. Consider a degenerate Fermi gas in an ax- 
ially symmetric trap slowly rotating around the symme- 
try axis. The density n s of the superfluid component can 
be defined as the rigidity of the superfluid with respect 
to an infinitesimal twist in the boundary conditions [14|, 
ipa(r) — * e~ m ^Vo-(r), where is the azimuthal angle 
with respect to the symmetry axis, Vv(r) = J2j ^ja- l Pj( r )i 
and ¥>j(r) are the single particle wave functions. The 
twist generates a term — aj in the Hamiltonian, where 
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J is the current operator. The resulting supercurrent 
J s = n s a/m can be expressed in terms of correlation 
functions (|13I14|1 using the standard linear response the- 
ory. We find n s = n/2, where n is the particle density 
in the normal state. A similar calculation for nonzero 
yields n s = n - the superfluid density in the gapped 
steady state is the same as in the BCS ground state. 
Since n s is the second derivative of the free energy with 
respect to a, the jump from n to n/2 at the critical point 
agrees with the second order character of the transition. 
The reduction in n s in the state with zero A™ is consis- 
tent with vanishing of the gap in this state [15J. 

We see that the superfluid density is finite in the gap- 
less steady state. This reflects the existence of pair cor- 
relations (|14|l between fermions, a property that in equi- 
librium has been directly linked to nonzero n s [lfij . This 
situation is similar to the phenomenon of gapless super- 
conductivity in metals [13, where the superconducting 
state is also characterized by non-vanishing anomalous 
Green's functions and zero spectral gap. A gapless state 
in superconductors is usually a consequence of strong per- 
turbations - as the perturbation is increased, the super- 
conducting metal first goes into a gapless state with finite 
n s before it becomes normal. 

A crude qualitative understanding of the dynamical 
transition to the gapless steady state can be derived 
from the following thermodynamic argument. After 
the change of coupling the initial state Q has energy 
£(Af, Ai) > relative to the ground state with gap Af. 
In thermal equilibrium the system would have the same 
energy at a certain temperature To (A/, A^). Let us keep 
Ai fixed and vary A/. At Af = Ai, To — 0, while the 
critical temperature T c w 0.57A/ > 0. As Ay increases 
both To and T c grow, but T c grows faster and To never 
catches up with it. On the other hand, when Af de- 
creases, T c also decreases, while To grows. Evaluating 
To (A/, Ai), to the first order in the coupling constant we 
find that T > T c for A// A, < 0.52 + 0.17A and T < T c 
otherwise. Here A = g/d is the dimensionless BCS cou- 
pling constant and d = (ej+i — €j) is the mean level 
spacing. We see that decreasing the coupling beyond a 
certain value provides enough energy for the transition 
to the normal state to occur in thermal equilibrium. 

Let us discuss experimental manifestations of the dy- 
namical transition to the gapless steady state. After the 
initial sweep (gi — > <?/) on the BCS side of the resonance 
such that Af/Ai < e _7r / 2 , the gapless state is reached 
on a timescale ta = l/Aj. The pairing gap measured 
in Ref. Q is the energ y co st AT of breaking a Cooper 
pair. This corresponds [12] to removing the spin closest 
to the Fermi level in The minimum energy cost is 
AT = min|b(e) • s(e)| = A oo cos6'(0). Using Eq. ©, we 

obtain AT = A^^Jl- 41n 2 (A i /A / )/7r 2 . In the steady 
state, AT = 0, i.e. no gap in the rf absorbtion spectrum 
will be observed. Molecular condensate fraction 111 |2| 



vanishes in the gapless state, since Nk—o/N cx A 2 ^ [lij . 

In conclusion, we determined the dynamics of the 
paired state of cold atomic fermions following an abrupt 
lowering of the pairing strength, </j — * gf. On a short 
1/Ai timescale, where A^ is the equilibrium gap for the 
initial coupling gi, the system goes to a non-stationary 
steady state with a constant gap A^ and a distribution 
function given by Eqs. I|8I9|) . When the coupling is re- 
duced so that Af/Ai is smaller than the critical value of 
e _7r / 2 , the steady state gap vanishes. The decay law of 
the time-dependent order parameter A(t) changes from 
power law to exponential as A f / Ai goes through the crit- 
ical point, Eqs. (|11I12|I . The gapless state combines fea- 
tures of normal and superfluid states. In particular, the 
gap and the condensate fraction vanish, while the super- 
fluid density is nonzero. 

We thank E. Abrahams, B. Altshuler, and J. 
Schmalian for useful discussions. M. D. was in part sup- 
ported by the U.S. Dept. of Energy, Office of Science, 
under Contract No.W-7405-ENG-82 and by the ICAM. 
After this work was completed, we became aware of a 
work by Barankov and Levitov|l9j that reports some re- 
sults similar to ours. 



[1] C. A. Regal, M. Greiner and D. S. Jin, Phys. Rev. Lett. 

92, 040403 (2004). 
[2] M. W. Zwierlein et.al, Phys. Rev. Lett. 92, 120403 

(2004) . 

[3] J. Kinast et. al. Phys. Rev. Lett. 92, 150402 (2004). 

[4] C. Chin et.al., Science 305, 1128 (2004). 

[5] R. A. Barankov, L. S. Levitov, and B. Z. Spivak, Phys. 

Rev. Lett. 93, 160401 (2004). 
[6] E. A. Yuzbashyan, B. L. Altshuler, V. B. Kuznetsov, V. 

Z. Enolskii, J. Phys. A 38, 7831,(2005); Phys. Rev. B 72, 

220503(R) (2005). 
[7] M. H. Szymanska, B. D. Simons, and K. Burnett: Phys. 

Rev. Lett. 94, 170402 (2005). 
[8] G.L. Warner and A. J. Leggett, Phys. Rev. B 71, 134514 

(2005) . 

[9] E. A. Yuzbashyan, O. Tsyplyatyev, B. L. Altshuler, Phys. 

Rev. Lett. 96, 097005 (2006). 
[10] K. M. O'Hara et. ai, Science 298, 2179 (2002). 
[11] T. Loftus, C. A. Regal, C. Ticknor, J. L. Bohn, D. S. Jin, 

Phys. Rev. Lett. 88, 173201 (2002). 
[12] P. W. Anderson, Phys. Rev. 112, 1900 (1958). 
[13] A. F. Volkov and Sh. M. Kogan, Sov. Phys. JETP 38, 

1018 (1974). 

[14] M. E. Fisher, M. N. Barber, and D. Jasnow, Phys. Rev. 
A 8, 1111 (1973). 

[15] J. R. Schrieffer, Theory of superconductivity, W. A. Ben- 
jamin, New York, 1964, Ch. 8. 

[16] W. Kohn and D. Sherrington, Rev. Mod. Phys, 42, 1 
(1970). 

[17] K. Maki, in Superconductivity (edited by R. D. Parks), 
Marcel Dekker, New York, 1 969, v.2, p 1035. 

[18] R. B. Diener and T. L. Ho, cond-mat/0404517 E. Alt- 
man and A. Vishwanath, Phys. Rev. Lett. 95, 110404 



(2005). 

[19] R. A. Barankov and L. S. Levitov, cond-mat/0603317 



